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Abstract 

The Picard-Fuchs equations for N = 2 supersymmetric SU(N C ) 
Yang-Mills theories with massless hypermultiplets are obtained for 
N c = 2 and 3. For SU{2) we derive the non-linear differential equa- 
tions for the prepotentials and calculate full non-perturbative correc- 
tions to the effective gauge coupling constant in the weak and strong 
coupling regions. 
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1 Introduction 



The work of Seiberg and Witten |], || on the exact solution to the low-energy 
effective action of N = 2 supersymmetric Yang-Mills theory has afforded not 
only a renewed insight into the charge confinement |J and the chiral sym- 
metry breaking ||, but also a marvelous insight into topological invariants 
[[| and conformal field theories in four dimensions M. The key ingredients 
in obtaining these exact results are duality and the appearance of mass- 
less monopoles/dyons in the strong coupling regions of the theory. In the 
weak coupling region, the exact solution enables us to determine full non- 
peturbative instanton corrections to the effective coupling constant, whose 
evaluation is otherwise quite cumbersome in the standard framework of quan- 
tum field theory. 

The low-energy effective action is described in terms of a prepotential 
which is a single holomorphic function of superfields of N = 2 £7(1) vector 
multiplets. The exact solution for the prepotential may be characterized by 
the period integrals of the special type of one-form on a hyper-elliptic curve. 
The curves associated with a variety of iV = 2 supersymmetric Yang-Mills 
theories and QCD have been studied extensively [j5| — pOp. The moduli space 
of the curve contains singularities, at which some solitons become massless. 
In order to investigate the strong coupling physics, one needs to evaluate the 
period integrals near the singular locus. An efficient approach to study this 
problem is to use the Picard-Fuchs equation, the differential equation which 



the periods obey [|TT], y, [lj, [14 



In the present article, we shall investigate the quantum moduli space of 
N = 2 supersymmetric Yang-Mills theories with massless hypermultiplets 
and gauge groups SU (2) and SU (3). We shall derive the Picard-Fuchs equa- 
tion for the scalar part of N = 2 U(l) vector multiplets and their duals. By 
solving the non-linear differential equation obeyed by the prepotential for 
G = SU(2), we will explicitly evaluate the non-perturvative contributions in 
the prepotential in both weak and strong coupling regions. 
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2 N = 2 Supersymmetric Yang-Mills Theory 
with Massless Hypermultiplets 

We begin with reviewing some basic properties of the low-energy effective 
action of the N = 2 supersymmetric G = SU(N C ) QCD 0,0. In JV = 1 
superfield formulation, the theory contains chiral multiplets $ a and chiral 
field strength W a (a = 1, . . . , dimG) both in the adjoint representation of 
G, and chiral superfields Q l in the N c and Q l (i = 1, • • • , Nf) in the N c 
representation of the gauge group. The superpotential is given by W = 
QT a <& a Q J r M l jQ i Qi where T a is the generator of G and Mj is the mass matrix. 
Along the flat direction the scalar fields of $ get vacuum expectation values, 
which break the gauge group to the Cartan subgroup U(l) r where r = N c — 1 
is the rank of G = SU(N C ). When the squark fields do not have vacuum 
expectation values, the low-energy effective theory is in the Coulomb branch 
and contains r U(l) vector multiplets (A 1 , W^) (i = 1, ■ • • , r), where A 1 are 
N — 1 chiral superfields and W % a are N = 1 vector superfields. The quantum 
moduli space may be characterized by the low-energy effective Lagrangian L 
with the prepotential T{A) 

L = — Im (J d 2 6d 2 6A Di A l + d^r^W^W^ , (1) 

where Am = J^j is a field dual to A 1 and = d %Q Ai the effective coupling 
constants. We denote the scalar component of A\ A Di by (3*, QiDi- The 
pairs (ar>i,o>i) are the Sp(2r, Z) section over the space of gauge invariant 
parameters Sj (i — 2, • • • , N c ) defined by det(x — 0) = x Nc — Ya=2 SiX Nc ~ % . 
The quantum moduli space of the Coulomb branch is parametrized by the 
gauge invariants Sj and the eigenvalues mi, . . . , m^. of the mass matrix. 

The sections (az>i,CLi) are obtained as the period integrals of a meromor- 
phic differential A over the hyper-elliptic curve C 0: 

V 2 = C(x) 2 -G(x), 

N c \2Nc-Nf Nf—Nc 

C(x) = x N°-Y, 8iX N°-* + - J2 U{m)x N ^-\ 

i=2 * i=0 



I 

i=l 



G(x) = A^-^n^ + m,), (2) 
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where A stands for the QCD scale parameter and ti(m) is defined by Y\i ! Ji(x+ 
m i) — YlfJoU{nT-)x Nf ~' 1 with to{m) = 1. The terms proportional to A 2Nc ~ N f 
in C(x) are absent in the case of Nf < N c . Notice that for N c = 2, Nf = 2 
theory C(x) should read 

C(x) =x 2 -s 2 + ^. (3) 
The meromorphic differential A is given by 

A = ^dtog^. (4) 
2m C + y 

This differential satisfies the relations 

— = -u Nc+1 -i + d{*), (5) 

where uoi = — X y dx is the basis of the holomorphic one-forms of the hyper- 
elliptic curve C. Choosing the basis of homology cycles {A, h B,j) on C with 
the canonical intersection form Ai ■ Bj = 5ij, Ai ■ Aj = Bi ■ Bj = 0, one can 
express the pairs (a,Di,cn) as 

Oii- a D = / A, 

JBi 

(Aj-Aj_i)-a = / A, i = l,---,r, (6) 

J Ai 

where are the simple roots of G and Aj the fundamental weights of G 
(A = 0). 

The BPS mass formula M > \qia,i + hiar>i\ shows that the massless soliton 
may be characterized by the vanishing cycle + hiBi ~ 0. Therefore the 
singularities of the quantum moduli space are determined by the discriminant 
A of the curve. Let us enumerate explicit forms of the discriminants for 
the SU(N C ) theories (N c = 2, 3) with N f (< 2N C ) matter hypermultiplets. In 
the following denotes the scale parameter corresponding to the iV^-flavor 
theory. 

• N c = 2 (u = s 2 ) 

A 2 = 256Aq(-u - Aq)(m + Aq) 
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A 2 = -A^(27A? + 256m 3 ) 

A 2 = A 4 (8«- A 2 ) 2 (8« + A 2 ) 2 /16 

A 2 = A 2 3 (A 2 3 - 256w)w 4 

(7) 

• iV c = 3 (u = s 2 , v = S3) 

Aq = 64A^ 8 (-27(t; - A 3 ) 2 + 4w 3 )(-4w 3 + 27(v + A 3 ) 2 ) 
A 3 = (-3125AJ 5 + 256A^m 5 + 22500Aj°uu - 102Au% 

-43200AfuV + 13824uV - 46656w 5 ) 
A 3 = 64A^V(-4(A 2 - M ) 3 -27t; 2 )(-4(A 2 + M ) 3 + 27t; 2 ) 
A 3 = A^(-A 3 + 4v) 3 (729A^ 3 - 256M 6 - 3888A 3 M 3 t; 

+3456wV + 2916A 3 u 3 - 11664w 4 ))/16 
A 3 = 4A^ 4 (A 4 m - 8A 2 w 2 + 16m 3 + 2A 3 w - 72A 4 uv - I08v 2 ) 

{-A\u + 8A 2 m 2 - 16m 3 + 2A 3 w - 72A 4 wv + 108w 2 ) 
A 3 = 4A 3 v 5 (A 2 m 5 - 256w 6 + 2A 3 w\; - 528A 5 tA; + A\uv 2 

-300A 2 wV + 3456uV + 4A 3 w 3 - 4212A 5 w 3 - 11664v 4 ). 

(8) 

Note that for the gauge group 577(2) the zeros of the discriminant realize 
the discrete Z^N f symmetry which arises from the U(1)k anomaly and the 
parity transformation in the flavor symmetry 0(2Nf). Order of the zeros 
corresponds to the number of massless hypermultiplets appearing at the sin- 
gularities. Massless monopoles obey the spinor representation of the flavor 
symmetry SO(2Nf). For SU(3) symmetries are less obvious particularly for 
odd flavors. 

3 The Picard-Fuchs equations 

Let us study the Picard-Fuchs equations for (am, a») and examine their solu- 
tions. We first discuss the SU (2) case in detail for simplicity. One can change 
the curve of fourth order to the third one by using the Mobius transformation. 
The curves are given by 

y 2 = X 2 (X-U) + ^AqX (9) 
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for Nf = and 



y 2 = x 2 (x-u)-±A^- N f\x-u) N f-i 



(10) 



for Nf = 1, 2, 3. The meromorphic differential A for this curve is given by 

v/2 2u - (4 - N f )x 



X 



dx 



in 



(11) 



which satisfies 



dX \pldx ,. . 

— = v dM. 

du 8tt y 



(12) 



It is not difficult to see that the period EE = § X satisfies the differential 
equation 

f ,d 2 U 

p(«)^ + n = o, (is) 



where [11, 12, 14 



p(u) 



du 2 



f 4(« 2 - A*) 
„ 2 27A? 

2 A 2 



N, = 



4 m 



A 6 ! 



N 



f 



I V 64/ 



(14) 



Note that the regular singularities of (p~3|) correspond to the zeros of the 
discriminant. Introducing a new variable 



256 / it 



64 ^ 

256 



27 I Af 

/ \ 2 ' 
I U \ 



hi) 

A| 



JV> = 
iV> = l 
JV> = 2 



(15) 



N 



f 
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it is shown that the Picard-Fuchs equation ( p~3| ) turns out to be the Gauss's 
hypergeometric system 

d 2 Il dTl 

z(l - z)-^ + (7 - (a + (3 + - a(3U = 0, (16) 

where 

a = p = w^wr 7 = i^- (17) 

Using the fundamental solutions of fllB]) , we may evaluate the power series 
solutions for (ap, a) near the singularities, which has been explicitly done in 
our previous work |L4j]. 

Next we write down the Picard-Fuchs equations for SU(3) theory with 
flavors Nf < 5. The periods II = § A now depend on two parameters u and 
v . The Picard-Fuchs equations take the form 

£ill = 0, £ 2 n = o, (18) 

where C\ and £2 are given by 

• N f = 

£1 = -(-27A^ + 4m 3 + 27v 2 )d 2 + 12uvd u d v + 3vd v + 1 
3 

£ 2 = -(-27A{j + 4u 3 + 27v 2 )dl + 12uvd u d v + 3vd v + 1 (19) 
u 

• N f = l 

-25A?m 2 + 84u 3 v + 405m o9 64m 4 - 1125A?v + 2160m; 2 n n 

Ci = dl H d u d v 

45w v 180v 

16m 3 + 135v 2 , 

+ 45, + 1 

4(-25Afw 2 + 8Au 3 v + A05v 3 ) 2 
-25Af + 84m; 

(-25Af + 36w)(-25Af + 96w) 
+ -25Af + 84m " 4 

3w(-25Af + 36m) 



-25Af + 84m 



d v + 1 (20) 
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-8A 4 m + 8m 3 + 27v 2 2 4(2A» - 4A 4 m 2 + 2m 4 + 27uv 2 ) 

3 " + 9^ u v 

-8A 4 w + 8m 3 + 27v 2 . 

a„ + 1 



9w 

-8 A 4 m + 8u 3 + 27v 2 2 3v (A 4 + 3m 2 



OS + — ^ + 1 (21) 



2m u u 



{4v- A 3 )(4m 3 + 9^ 2 ) a2 m(64m 3 - 81A 3 t; + 432t; 2 ) 
4^ ^ + 36^ 

16m 3 + 27w 2 

+ 9, * + 1 

4^ M (-9AI + 32,) aA _^ + 1 

M 4 



-Aju - 56A 2 w 2 + 240m 3 + 324v 2 a2 -8A 2 m 2 + 32m 3 + 27v 2 

36 ^ + 9, ^ 

2AV - 16A 2 m 3 + 32m 4 - 9AV + 108mw 2 

H „ + 1 

9f 

-A 4 m - 56A 2 m 2 + 240m 3 + 324v 2 2 
A 2 + 60m u 
^(A 2 + 12m)(-A 2 + 36m) 3t;(A 2 -36M) 
+ A 2 + 60m + A 2 + 60m + 



(23) 



w(-2A 2 m 2 + 528m 3 + 57A 5 mt; + 324w 2 ) g2 
5A 5 m + 36v " 
-5A 2 m 3 + 1280m 4 - A|mm + 332A 5 M 2 t; - 4A 2 t; 2 + 1728m?; 2 
+ 4(5A 5 m + 36v) 



-5A 2 u 2 + 1280u 3 + 88A 5 w + 432w 2 „ 

-I 8 4- 1 

4(5A 5 w + 36w) v + 

-2A 2 w 2 + 528m 3 + 57A 5 uv + 324v 2 o2 

C - = — 132i * 

-A 5 w 3 - 5A?uv + 8QAu 2 v - 36A 5 v 2 n n A 5 w + 324w n 

-9 u o v — d v + 1. 



132m 132 

(24) 

The Picard-Fuchs equations for N c = 3, Nf = theory have been investigated 
in detail by Klemm, Lerche and Theisen fl2"|| . They show that the system 
reduces to the Appell's hypergeometric system of type F4. Although other 
cases are not classified as the Appell's system, the solutions of these equations 
may be obtained in the form of the power series expansion in u and v. 



4 Calculation of the Prepotentials 

Solving the system of the Picard-Fuchs equations, we evaluate the non- 
perturbative corrections to the prepotentials, which can be regarded as the 
multi-instanton contributions in the weak coupling region and the threshold 
corrections in the strong coupling regions. In the previous work [TJ[] , we have 
calculated some non-trivial corrections by using the solutions of the hyper- 
geometric system. In this article, we shall instead take more direct approach 
by solving the non-linear differential equation for the prepotential in SU(2) 
theory. The present approach was proposed originally by Matone for 377(2) 
pure Yang- Mills theory |13| . 

We start with the Picard-Fuchs equation 



13). Notice that the equation 



carries no first derivative term. This implies that the Wronskian of the system 



W 



a ajj 
d u a d u aD 



(25) 



is independent of u; d u W = 0. When we integrate this with respect to u there 
appears a constant, the explicit value of which is evaluated at any singularity 
in the moduli space. One finds that 



W 



ib 
4V 



(26) 
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where 6 = 4 — Nf is the coefficient of the beta-function of 377(2) N = 2 
QCD. In the present case, we can further simplify the equation ( p6|) since 
the Wronskian is written in the form of the total derivative 



W = d u (a^ - 2T 



(27) 



Thus we get the following relation between the prepotential and the gauge 
invariant parameter u 



a— 27 

oa 



ib 
An 



u + const. 



(28) 



The constant term may be absorbed in the definition of T by shifts, which 
does not affect the form of the low-energy effective action. This identity has 
been generalized for any N = 2 Yang-Mills theories with or without matter 
hypermultiplets PT], and is important in view of its relation to the soliton 
theory |23| 



The relation (|28l ) allows us to express u in terms of a 
u = Q{a)^(a d -^-2A. 



ib \ da 



(29) 



The Picard-Fuchs equation becomes 

d 2 g 



-v{Q{a)) 



do? 



da J 



0. 



(30) 



This non-linear equation is utilized to determine the instanton coefficients 
recursively. Near the singularity at u = oo, the prepotential takes the form 



t |21og^ + Er=o^(0) (^) 
t I ( 4 - N f ) In 4- + IXo MNj ) ( ^ 



N f = 



J1 JV f 



N f 



1,2,3. 



(31) 



The first four coefficients are listed in the table 
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N f 


Fx 


^2 


^3 


^4 





i 

2 5 


5 

2 14 


3 

2 18 


1469 

2 18 


1 


3 
2 12 


3-17 

" 2 288 


5-7-11 

3-2 38 


3 2 -5-7-9679 

2 59 


2 


1 

2 11 


5 


3 
2^ 


13113 

2 55 


3 


1 

95 


1 

2 24 





5 
251 



and summarized by the formulas (Nf = 1, 2, 3) [} 
1 1 



F 2 (N f ) 
F 3 (N f ) 
T A (N f ) 



~2 b c(N f )2b 2 ' 

1 5 - 86 + 46 2 
~2 2b c{N f ) 2 ' 

1 (1 - 6) (1 - 26) (-23 + 306 - 166 2 ) 
2 Sb c(N f ) 3 172866 ' 

1 (1 - 26) (-677 + 39106 - 81246 2 + 84566 s - 46726 4 + 11526 5 ) 



2 4b c(N f Y 



294912( 



(32) 



where c(N f ) = rj Nf 2 8 (4 - N f ) N f- 4 (Nf = 1, 2, 3) and Vl = -1, 7] 2 = r] 3 = 1. 
The Nf = result [|l|, |13| may be obtained by replacing c(Nf) with 1 in the 
Nf = 2 case. The results ([!2|) indeed agree with those obtained previously 

Similarly we may solve the non-liner differential equation in the strong 
coupling region where monopoles become massless. In this case the dual field 
an is a good coordinate around the singularity u = uq where 



u 



A 2 


N f 


= 


-3-2~§A 2 


Nf 


= 1 


A 2 /8 


Nf 


= 2 





N f 


= 3 



(33) 



*We have corrected an error of ^(Nf) in the previous work 
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Introduce the dual prepotential J-d{o>d) satisfying a = ^s.. Then, instead 
of (p7D , the Wronskian is expressed in terms of To- Hence around u = uo we 
get 

Similarly to the weak-coupling computation we now obtain the non-linear dif- 
ferential equation which makes it possible to calculate Tb recursively. Using 
the expansion 

T D = ±-al (Hog ($-\ + £ F Dn {N f ) (%>X\ (35) 



with k — I (Nf — 0) and fc = 2 7V f~ 1 (JV> = 1, 2, 3), we find the coefficients as 
follows 

^b_i(JV» = l6kb 2 c(N f ), 

2 -2k (l-26)(-5 + 26) 



2 (iV/) 
^ s(JV» 



3c(iV / ) 166 2 

-2k (1 - 26)(-61 + 1906 - 92& 2 + 8& 3 ) 
2c(JV/) 2 11526 4 ' 

2 • 2k (1 - 26) (1379 - 81626 + 145966 2 - 72886 3 + 9606 4 ) 
5c(iV / ) 3 1105 9 26 6 ' 

(36) 



where c(l) = -3" 1 / 2 • 2~ 17 / 6 ,g(2) = -i2~ 9 / 2 and c(3) = 2" 13 / 2 . The Nf = 
result is obtained by putting 6 = 2 and c(0) = — z/4. Here we have corrected 
some misprints in the previous work W% . 



For SU(3) theory the prepotential T obeys the basic identity analogous 
to ( p8l) [pi] , p2| . In contrast to the SU(2) case, combining this with the 
SU(3) Picard-Fuchs equations does not yield a simple recursion relation. It 
is not clear at present if the existence of such recursion relation is peculiar to 
SU (2) theory. It will be interesting if one could find an additional constraint 
equation for T which eventually leads to the direct evaluation of instanton 
corrections. 
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